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1 Introduction 



One of the long standing problems in multiparticle dynamics is the descrip- 
tion of multiplicity distributions and correlation functions within a common 
formalism which can emphasize the dynamical processes underlying multipar- 
ticle production. In this talk, I will illustrate one step in this direction, by 
discussing how the sign oscillations of the ratio of factorial cumulant moments 
to factorial moments of the Multiplicity Distribution (MD) can be related to 
the dynamics of the early stage of partonic evolution in e^e~ annihilation. 

In Section ^ the relevant observables are defined and illustrated with ex- 
amples; in Section ^ and | theoretical calculations and experimental results 
are reviewed; in Section E the shoulder structure in the MD and the oscilla- 
tions of moments are related to hard gluon radiation; finally in Section ^ the 
effect of flavour quantum numbers on the moments of the MD in 2-jet events 
is discussed; conclusions are drawn at the end. 



2 Definitions and examples: the effect of truncation 

In this section the definitions of some differential and integral observables and 
their relationships will be recalled; see referenceEl for further details. The defi- 
nitions, summarized in Table |l|, concern n-particle distributions: the exclusive 
ones, P„(yi, . . . , ?/„) where all n particles produced at rapidities yi, ■ ■ ■ , j/n are 
observed, and the inclusive ones, Qn{yi, ■ ■ ■ ,yn) where at least n particles are 
observed. These quantities are of course related: 

Pn{yi, ■■■ ,yn) ^ Qn{yi, ■ ■ ■ ,yn) + 

j— / dy[...dy'„^Qn+m.{yi,--- ,yn,y'i,--- ,y'„,) (1) 

^ — ' TO / 

m=l 

From the inclusive distributions, via cluster expansion, one gets the correlation 
functions C„(?/i, . . . ,?;«), thus subtracting from Qn the statistical, uninterest- 
ing correlations due to combinations of lower order distributions Qn-i ■ ■ - Qi] 
howeverJ;here remain correlations related to fiuctuations in the number of par- 
ticles nn it is therefore only in an approximate sense that one can say that 
when C2 is positive particles like to cluster together (which one relates to the 
dynamics of the process) and when C2 is negative particles like to stay away 
from each other (which one relates to the effect of conservation laws) . 

By integrating the differential observables just discussed one obtains the 
integral observables which will be the subject of the rest of this talk (see 
Table 0): the multiplicity distribution, the factorial moments and the factorial 
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Table 1: Definitions of relevant differential (on the left) and integral (on the right) quantities; 
a is the inclusive cross section, (Tinel the inelastic cross section, yi is the rapidity of the i-th 
particle. The integrals in the right-hand column are over the full phase space. 



Exclusive distribution 


Multiplicity distribution 


Prr{V\,--- ,yn) 


-Pn = ^ y dyx... dynP„{yi,. .. ,y„) 


Inclusive distribution 


Factorial moments 


Qn{yi, ■■■ ,yn)- , , 

o-inci dyi ... dyn 


-Fn = J dyi. ..dynQ„(yi, ... ,?/„) 


Correlation function, e.g. 


Factorial cumulant moments 


C2{yi,y2) = Q2{yi,y2) - Q\(y\)Q\{y2) 


K„ = J dyi . . .dy„C„{yi, . . . ,y„) 



cumulant moments. They are linked by the foUowing relationships, analogous 
to Eq. 1, which allow to obtain the moments directly from the MD: 

C30 

^^« = ^r-(r-l)---(r-n+l)P, (2) 

r—n 

( — l\ 

Kn=F^-Y,r jK^^rFr (3) 

It is apparent from the definitions that Fn and Kn receive contributions from 
events with at least n particles, which means that moments of high order are 
very sensitive to the tail of the MD. In particular, as will be shown in the 
following, it is interesting in this respect to study the behaviour of the ratio 

H,^§^ (4) 

as a function of the order q: it is qualitatively different for different distribu- 
tions and turns out to be suitable for analytical calculations (see Section H). In 
order to illustrate these properties. Figure [l] (left column, dashed lines) shows 
the ratio Hg for some of the most common discrete distributions. It is worth 
pointing out that for the Poisson distribution, which one usually associates 
with a lack of dynamical correlations, the ratio Hq is zero for g > 1; for the 
geometric and for the negative binomial distributions, which one usually as- 
sociates with dynamical effects because they give rise to positive correlations, 
the ratio Hq is always positive, but decreases toward zero; for the binomial 
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Figure 1: The ratio Hq for the most common discrete distributions. All the distributions are 
chosen to have the same average number of particles (namely, n = 10); other parameters, if 
any, are listed in the figure. The plots in the left column correspond to the full distributions 
(dashed lines), and to the even-only distributions (solid lines joining the diamonds). The 
plots in the right column correspond to the same distribution after truncation has been taken 
into ppcount as in Eq. ^; the effect of the truncation changes with the highest multiplicity 
ricutfl the values chosen here are such that the discarded part is always less than 1% of the 

cumulative distribution function. 
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distribution, which one usually associates with conservation laws because of 
negative correlations, the ratio Hq changes sign according to the parity of q. 
Before applying these considerations to the data in full phase space, one should 
remember that only even multiplicities are allowed for charged particles. In or- 
der to show how the suppression of the odd multiplicities affects the moments, 
the left column in Figure |l| presents also the ratio Hg vs the order q subject 
to the condition that P„ = if n is odd (solid lines and points): it is seen that 
there is a small distortion with respect to the values of the full distributions. 

As already mentioned, the ratio Hq probes the tail of the MD. Unfortu- 
nately, the tail is the most difficult part to measure experimentally, because 
only a finite number of events can be collected. This results in a MD which is 
truncated at some point: 

P„ocj^" (5) 
I otherwise 

where the proportionality factor ensures proper normalization. The moments 
of high order are affected: the factorial moments are smaller than in the full 
distribution, and one finds -that the factorial cumulant moments oscillate in 
sign as the order increasesS An example of this effect can be seen again in 
Figure]^; its importance for experimental MD's will be discussed in Section |[ 



3 The ratio Hq in perturbative QCD 

The first suggestion that the ratio Hq could be usefuLiiJ the study of MD's 
came from analytical calculations in perturbative QCDpQ where it was shown 
that one expects oscillations of the ratio Hq as a function of the order q. 
In the followmg I will summarize the derivation in the framework of pure 
gluodynamicsp starting from the evolution equations for the generating func- 

[g{y + \nx; z)giy + ln(l - x); z) - g{y; z)] (6) 

where the first term in parentheses is the DGLAP kernel for the process 
g gg with the emitted gluon carrying a fraction x of the parent momentum, 
y = ln(fcj^/(5o) is the evolution variable {k± the jet's transverse momentum 
and Qq the cutoff, which limit the integration interval via x{l — x)kj_ > Qo), 
(^s{y) is the running coupling constant and is the number of colours. No- 
tice that the terms in the square brackets take recoil into account, thus going 
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beyond the pure Double Log Approximation (DLA). In order to solve Eq. 6 for 
the moments of the distribution one expands it in powers of 3^, remembering 
that, when calculating correlations of order q, one should use qj as expansion 
parameter (because the NLO contribution to such correlations is proportional 
to 97, rather than 7; 7 is here the anomalous dimension). Assuming asymp- 
totic KNO scaling, so that only the average number of partons depends on 3^, 
while the normalized moments of higher order are constant as 3^ — > cx), one 
calculates now the derivatives in z = of the generating function Q, obtain- 
ing the factorial moments; the factorial cumulant moments are then computed 
via Eq. 3. The normalization of these moments is not fixed but can be elimi- 
nated by taking their ratio. One finds for the ratio Hq a negative minimum at 
q 5 and sign oscillations at larger q. This result has been qualitatively im- 
proved hy including an estimate of the contribution of the vertices q qg and 
g qq-p further confirmation comes from the exact solution of the full evolru 
tion equations including both quarks and gluons in the case of fixed couplingQ 
It could be noted that the oscillations found in this last case happen around 
the (monotonically decreasing) value of Hq of a negative binomial distribution 
(NBD) with A: « 5. On the other hand, it is interesting to point out that in 
DLA the oscillations disappear and the ratio Hq is very close to that of a NBD 
with fc « 3. 

However, before doing comparison of partonic results with experimental 
data, one should, investigate what the possible role of hadronization is. Here 
I will just recall^ that if the simplest possible Generalization to Local Parton 
Hadron Dualitylj (GLPHD), which requires the proportionality of all inclusive 
distributions at parton (p) and hadron (h) level via a single parameter p: 

Q^(2/i,... ,y„)=p"QP(yi,... ,y„), (7) 

is used, it is easy to show, by direct substitution into the formulae of Table 0, 
that one obtains: 

h:i = ^=h^ (8) 

Thus, if the GLPHD hadronization prescription is used, the parton level result 
on the ratio Hq can be directly applied to the hadrons: in this case, it is seen 
that the QCD prediction for e+e^ annihilationQ fails quantitatively. 

4 The ratio Hq in experiments 

The theoretical work described in the previous Section triggered the analysis 
(a posteriori) of available data on MD's in order to extract the features of the 
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ISR pv \fs = 30.4 GeV 

n = 10.7 k = 11.0 ncut = 26 



ISR pp 52.6 GeV 

n = 12.2 fc = 9.4 ricut = 32 



ISR pp^/s^ 62.2 GeV 
n = 13.6 k = 8.2 ricut = 38 



UA5 pp^= 200 GeV 
n = 21.2 fc = 4.8 Ticut = 58 



UA5 pp yfs = 546 GeV 

n = 28.3 fc = 3.7 n^ut = 100 



UA5 pp = 900 GeV 

n = 35.2 k = 3.7 ricut = 104 



Figure 2:[— iThe ratio Hq vs the order q for a collection of experiments in hadron-hadron 
collisions|i3 The solid line joins data points and is drawn to guide the eye only. The dashed 
line shows the prediction of the NBD which fits well the MD's, after the truncation effect 

has been taken into account. 



ratio It should be immediately said that this is not an easy task: in order 
to meaauxe the MD in a modern experiment an 'unfolding' procedure has to 
be usedlHl the resulting published MD has correlations between adjacent bins 
which are not taken into account when moments are extracted from it without 
knowledge of the correlation matrix. This makes the calculation of the errors 
on the moments thus obtained very difficult: it should be considered as an order 
of magnitude estimate rather than a precise determination. Keeping this in 
mind, one can proceed to review some experimental results. The main point 
here is that oscillations of very different amplitudes are seen in the data in all 
reactions, in most cases (but not all! see below) compatible with being due to 
the truncation of the MD. As an illustration, in Figure || some results relative 
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to hadron-hadron collisions arc shown: large oscillations arc found. The MD's 
from which these moments are extracted are well described (with the exception 
of UA5 data at yfs — 900 GeV) by a NBD. Also shown in the same figure is 
the ratio Hq predicted by the fitted NBD after taking into account the effect 
of truncation (as in Eq. ^) : it is seen that the oscillations thus obtained are of 
the same order of magnitude as the data, so that no further dynamical effects, 
bempnd those described by a single NBD, are apparent (see also the analysis 
in§). 

The case is different in e+e^ annihilation data at the ^^-jpeak, as exem- 
plified in Figure ^, where data from the SLD CoUaborationtil are compared 
with the predictions of a NBD (dotted line) and a truncated NBD (dot-dashed 
line). It should be mentioned that in this plot the errors on the data are the 
statistical errors as published by the SLD Collaboration, which take into full 
account all correlations from the unfolding matrix. Clearly the truncated dis- 
tribution cannot describe the data: it will be shown in the next section how 
one can relate these oscillations to hard gluon radiation. 



5 Common origin of the shoulder structure and of the oscillations 

A very interesting feature in the data on charged particles MD's at the peak 
is the shnuldep^tructure which is clearly visible in the intermediate multiplic- 
ity ranget§ll30 This peculiar behaviour appears evident if one looks at the 
residuals (difference between the data and fitted values, in units of standard 
deviations) with respect to a NBD (or even a Lognormal Distribution) liJ the 
curve starts below the data, goes above and then below again. A pattern is 
seen, instead of a random sign and size of xesiduals. 

The DELPHI Collaboration has shownta that the shoulder structure in the 
MD in e^e~ annihilation can be explained by the superposition of the MD's 
coming from events with 2, 3 and 4 jets, as identified by a suitable jet-finding 
algorithm, and that the MD's in these classes of events are well reproduced 
by a single NBD for a range of values of the jet-finder parameter, j/min- It is 
thus suggested that the shoulder is associated with the radiation of hard gluons 
resulting in the appearance of one of more extra jets in the hadronic final state. 
One should also recall that a shoulder structure similar to the one obseipasd in 
e^e~ annihilation has been observed in pp collisions at high energiesllj'c3 and 
was shown to be well described by a S-patameter parametrization in terms of 
the weighted superposition of two NBD'sEl 

Following these observations, a parametrization of the MD which is the 
weighted sum of two components, one to be associated with S-jei. events and 
one to be associated with events with 3 or more jets, was proposedE3 The weight 
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Figure 3: The ratio of factorial cumulant moments over factorial pwjments, Hq as a function 
of q; experimental data (diamonds) from the SLD CoUaborationLj are compared with the 
predictions of several parameterizations, with parameters fitted to the data on MD's: a full 
NBD (dotted line); a truncated NBD (dot-dashed line); sum of two full NBD's as per Eq. M 
(dashed line); sum of two truncated NBD's as per Eq. nd (solid line). 



in this superposition is then the fraction of 2-jet events, which is experimentally 
determined, not a fitted parameter. This decomposition depends of course on 
the definition of jet, and in particular on the value of the parameter which 
controls the jet-finding algorithm. The DELPHI CoUaborationEa has used the 
JADE algorithm and published values for the 2-jet fraction and for the MD's 
at ymin = 0.02, 0.04, 0.06, 0.08. 

As for the particular form of the MD in the two components, the NBD 



was chosen because it has successfu 
of events with fixed number of jets 
MD's with a four parameter formula: 
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been fitted to the data for the samples 
In practice a fit was performed to the 



P aP^^^{ni,ki) + {l-a)P^''''{n2,h) if n is even 
I otherwise 

Here P^^^{n, k) is the standard NBD of parameters n and k; notice that the 
charge conservation law is taken into account. The proportionality factor is 
fixed by requiring the proper normalization for P„. 

Results of the fit to the data of four experiments can be summarized as 
follows: the per degree of freedom are equal to or smaller than 1, and values 
of the parameters are consistent between different experiments; they are also 
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consistent with those obtained by the Delphi Collaboration in fitting their 2-jet 
and 3-jet data separately with a NBD. 

Once the fits to the MD's have been done, one can compare the experi- 
mental data on iJg's with the values obtained from the fitted MD's by using a 
formula that takes into account the truncation effect too: 



where nmin and rimax are the minimum and maximum observed multiplicity, 
and a proportionality factor ensures proper normalization. A very good agree- 
ment with the data is obtained, as exemplified in Figure || (solid line). Notice 
that it is not possible to reproduce the behaviour of the ratio Hg without tak- 
ing into account the limits of the range of the available data: this can be seen 
again in Figure ^, where the dashed line corresponds to Eq. ^ and does not 
agree with the data. 

It can thus be concluded that the observed behavior of Hg^s results from 
the convolution of two different effects, a statistical one, i.e., the truncation of 
the tail due to the finite statistics of data samples, and a physical one, i.e., the 
superposition of two components. The two components can be related to 2- 
and 3-jet events, i.e., to the emission of hard gluon radiation in the early stages 
of the perturbative evolution. Notice that as the energy increases, the number 
of components also grows so that, asymptotically, the oscillations should be 
washed out, in agreement with the DLA expectations. 

6 Effects of flavour in 2-jet events 

A simple check of the picture described in the previous section consists in look- 
ing at the behaviour of thfi ratio Hg in 2-jet events, using the MD givenJay the 
DELPHI Collaborationil3 one does not expect oscillations. It is foundlllJ that 
oscillations are present, although their amplitude is one order of magnitude 
smaller than in the full sample. If one looks at the MD itself in 2-jet events, 
one discovers that the residuals, with respect to a NBD, show a structure sim- 
ilar to the one seen in the full sample; furthermore- the oscillations in the ratio 
Hg cannot be described by truncating one NBDE3 The hint that there could 
be a substructure comes once again from DELPHI data:tj they showed that 
the MD in a single hemisphere in a sample enriched in e+e~ — > bb events is 
identical in shape to the MD in a sample without flavour selection, except for 
a shift of 1 unit: 




(10) 



ib) ^ pialt) 



(11) 
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This effect could be due to weak decays of B-liadrons Tliis result hints at 
a possible substructure in terms of events with heavy quarks and events with 
light quarks, each sample contributing to ±he MD with one NBD. One can 
therefore try a parametrization of the formcJ 

Pnini,n,,k) = a,P,fB(n,,fc) + (1 - a,) P^"" (fii , k) (12) 

where now at is the fraction of bb events at the Z*^ peak, measured at LEP to 
be approximately 0.22. Notice that, as suggested by the DELPHI result, the 
parameter k is the same in the two NBD's: this is therefore a fit with three 
parameters. The parameters of the fits and the corresponding x^/NDF are 
given in Table |2| for different values of the jet resolution parameter j/min! the fit 
has been performed taking the charge conservation law into account similarly 
to Eq. H A really accurate description of experimental data is achieved. Notice 
that the best-fit value for the difference between the average multiplicities in 
the two samples, Sbi, also given in Table ^, is quite large. This difference 
grows with increasing j/min, i-e., with increasing contamination of 3-jet events. 
By looking at the residuals one concludes that the proposed parametrization 
can reproduce the experimental data on MD's very well, as no structure, is 
visibleES Furthermore this parametrization describes well also the ratio HqcB It 
is also remarkable that only the average number of particles depends on flavour 
quantum numbers, whereas the NBD parameter k is flavour- independent. 

One can thus conclude that the examination of the behaviour of the ratio 
Hq as a function of q has allowed to link the final hadronic level with the 
flavour composition of the event. 



Table 2: Parameters and Xn^DF of the fit to experimental data on 2-jet events MD's from 
the DELPHI CollaborationPi with the weighted superposition of two NBD's with the same 
parameter k (Eq. [l2| ); the weight used is the fraction of bb events. Results are shown for 
different values of the jet-finder parameter j/min- 





Vmin = 0.01 


Vmin = 0.02 


Vmin = 0.04 




16.81±0.21 


17.22±0.15 


17.98±0.15 


fib 


20.26±1.71 


21.96±1.57 


23.61±1.64 


k 


124±51 


145±53 


120±33 


xVndf 


17.4/16 


12.6/16 


27.5/20 


Sbi 


3.44±0.83 


4.6±0.5 


5.6±0.5 
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7 Conclusions 

The ratio of factorial cumulant moments to factorial moments, Hq, is a good 
observable for exploring substructures in hadronic final states. It is possible to 
calculate its behaviour in perturbative QCD and, after making allowance for 
the effect of finite statistics in the data, to extract dynamical information about 
the first stages of the perturbative evolution. In fact, hard gluon radiation 
explains the substructures observed in the full sample of events in e^e~ an- 
nihilation at the Z° peak; flavour dependent properties explain the additional 
substructures observed in the sample of 2-jet events. It is also remarkable that 
in this analysis the most elementary substructures are well described by nega- 
tive binomial distributions (down to 2-jet events of fixed flavour). A final point 
should be made, as the fits discussed in this talk have been performed on the 
published data, and therefore could not take into account the full correlations: 
the hope is that the interesting results thus obtained will spawn experimental 
work on the original data. 
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